Elastic and dynamic form factors of an atomic nucleus in the shell model 
with correction for the center-of-mass motion 1 



^3 ■ A.Yu. Korchm 2 and A.V. Shebeko 3 

O 

o , 

CN| , National Science Center "Kharkov Institute of Physics and Technology" 

61108 Kharkov, Ukraine 

03 
I—* 

I/"} | Abstract 

Analytical expressions for the elastic and dynamic form factors (FFs) are derived in the 
^-j- ' shell model (SM) with a potential well of finite depth. The consideration takes into account the 

motion of the target-nucleus center of mass (CM). Explanation is suggested for a simultaneous 
shrinking of the density and momentum distributions of nucleons in nuclei. The convenient 
■ working formulae are given to handle the expectation values of relevant multiplicative operators 



o 



X 



in case of the Is — lp shell nuclei. 



1. It is known that nuclear SM wave functions (WFs) do not possess the property of translational 
invariance (TI). Several methods in earlier and recent studies of nuclear structure (see, e.g., [1 —4]) 
have been proposed to transform any WF into one which is translationally invariant. Of great 
interest among these methods is the projection procedure considered in [4]. Along with other 
attractive features shown in [4], this procedure enables a comparatively simple evaluation of the 
corresponding CM correction to the purely shell quantities (see Ref. [5]). 

The approach developed in [5] is extended here to calculate the cross sections of the elastic 
and quasifree electron scattering on atomic nuclei with single-particle (s.p.) configurations more 
complex than the Is 4 one. In particular, we pay special attention to the physical interpretation of 
a simultaneous shrinking of the density and momentum distributions of nucleons in a nucleus due 
to the employed separation of its CM motion. 

2. By definition, the elastic FF in question is 

F(q) = (<$> intr | exp^n - R)\ | $ intr ), (1) 
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while the dynamic FF can be written as in [5], 

oo 

= ^ J exp(-iar)S(q,T)dT, (2) 

— CO 

^ 1 "" * 1 

S(q,r) = (® in tr I exp[«(pi - A~ P)qm~ r] | $ int r), 
a = uj + q 2 /2m, 

where u(q) is the energy (momentum) transfer, m the nucleon mass, 

A ^ A 

a=l a=l 

the CM coordinate (total momentum) operator of the nucleus composed of A nucleons so that 
r a (Pa) the coordinate (momentum) operator for nucleon number a, and §i n t r the intrinsic WF of 
the nuclear ground state (g.s.). 

Following [4], we take as | $j„ tr ) the vector 

I ®intr) = (R = | <£>>[<<£> | R = 0)(R = I $)]- 1 / 2 , (3) 

for a given trial (approximate) WF $ that may be nontranslationally invariant (nTI). Here, a round 
bracket, | ), is used to represent a vector in the space of the CM coordinate only. 

In the harmonic oscillator (HO) model, where the Slater determinant | $) is "pure" in the space 
of the CM coordinate (the Bethe-Rose-Elliot-Skyrme theorem [6,7]), one has 

F(q) = exp(^)F (g) (4) 

with 

F (q) = ($ | exp(zgri) | $) (4') 

and 

S(q,T)=eM b ^-)S (q,r) (5) 

with 

So(q,T) = ($ | exp(«p 1 gr/m) | $) (5') 

Here r is the oscillator parameter, b = poq/m, p = t-q 1 . 

Result (4) is widely used in applications starting from the work [8]. Note also that Eqs.(4) 
and (5) are valid within the HO model, being independent of any specific way to separate the CM 
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motion (e.g., the Ernst-Shakin-Thaler prescription (3) that is equivalent to the so-called "fixed CM 
approximation" in case of finite nuclei (bound systems) ). 

Now, before finding some analogs of relations (4) — (5) with an arbitrary WF $ (in particular, 
the Slater determinant constructed of the s.p. orbitals in a potential well of finite depth, say, the 
Woods-Saxon or Hartree-Fock ones), we would like to trace the CM corrections of the density and 
momentum distributions p(r) and i](p) within the HO model. In this connection, let us recall the 
general definitions for these quantities of primary concern: 

p(r) = (<S> int r | S(h -R-r)\ $ int r) = (2tt)- 3 J exp(-iyr)F(q)dq (6a) 

and 

■qip) = (<5> in tr | 6@! ~ P/A - p) | $ intr ) = 

= (2n)~ 3 (t / m) 3 J exp(—ipqT/m)S(q,T)dq (66) 
For the ls-lp shell nuclei we find in the HO model, 

Fo(Q) = (1 " ^g 2 r 2 )exp(-ig 2 r 2 ), 

5 (g,r) = (l-^6V)exp(-ioV) 

Substituting these expressions, respectively, into Eq.(4) and Eq. (5) we get with the help of 
formulae (6): 

p(r) = JL. + \^~) exp(-^), (7o.) 

^--^(^H-^Jexrt-g), (76) 



A- I A- I 



The intrinsic distributions without any CM correction are 

Po (r) ee ($ | 6(h - r) | $> = 7r- 3 /V - 3 (l + exp(~), 

,2 „2 



(8a) 



rjoip) ee ($ | o(A - p) | $) = vt- 3 / 2 p 3 (1 + 2(A j )P ) exp(-^), (86) 

By comparing Eqs.(7) and (8) one can see that the density and momentum distributions are subject 
to the equal changes. It does not seem to be surprising if we invoke the well known symmetry 
between the coordinate and momentum representations of the HO model Hamiltonian. One should 
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Figure 1: The variation of p o (r)/po(0) (solid) and of p(r)/p (0) (dashed) with r/r 

point out that, apart from the "most symmetrical" nucleus with A = 4, these CM corrections are 
not reduced to a simple renormalization of the oscillator parameter (e.g., the change r — > f ). 

The dependences p(r) and p(r) calculated by formulae (7) and (8) with A = 16 are depicted in 
Fig.l In this context, we show the ratio 



We see that an additional correlation between the nucleons, incorporated into (inherent in) 
the intrinsic density p(r), gives rise to their redistribution between the nuclear shells (from the 
ls-shell to the lp-one, while the center of the lp-distribution is shifted toward larger r- values). 
This rearrangement of the nuclear interior is accompanied by a decrease of the nuclear density in 
its peripheral region. It implies the corresponding increase of the probability to find the nucleons 
in the central part of the nucleus. Remind that both the density distributions (DDs) p(r) and 
po(r) are normalized to unity. We interpret these features of the intrinsic DD p(r) as its shrinking 




versus the ratio 
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in comparison with the DD po(r) , which embodies the spurious CM motion modes. The same 
interpretation can be applied to the momentum distribution (MD) 77 (p) vs f]o{p). 

Such a simultaneous shrinking of the DD and MD becomes more tractable if one evaluate the 
respective r.m.s. radii and momenta. By definition, one has 



It is readily seen that 



(r 2 ) = J r 2 p{f)df, (p 2 ) = J p 2 r/(p)dp 

(r 2 ) = J r 2 p (f)df= (r 2 ) + (000 | R | 000), 
( P 2 ) = J p 2 Vo (p)dp= (p 2 ) + A~ 2 (000 I P I 000), 



(9) 



(10a) 
(106) 



where the vector | 000) is used to describe the "zero" (lowest-energy) CM oscillations with respect 
to the origin of coordinates. A complementary smearing of po(r) and i] (p) compared respectively 
with p(r) and r](p) is due to the nonphysical motion mode. 

Further, let us consider the commutation relations for the relative coordinates f' a — f a — R and 
the canonically conjugate momenta p' a =p a — P/A (a — 1, . . . , A): 



a)k, 0' a )l] = ^ H (l ~ A' 1 ) (k, 1 = 1,2, 3) 



along with the original ones: 



[(f a )k, (Pa)l\ = tSta (k, 1 = 1,2, 3) 



The corresponding uncertainties meet the equations (see, for instance, [9], p. 54, and also Suppl. A 
to this translation ) : 



($ I (Af) 2 I $)($ I (Ap) 2 I $) = (r 2 ) (p 2 ) > - A , 



$ I (Af'Y I $)($ I (Af ') 2 I <*>) = (r 2 )*(p 2 )* > 



9 (A- 1 



12) 



;i2') 



4 V A 

for any state | $) normalized to unity, where we have introduced the expectation values (r 2 )^ = 
($ I (f^) 2 I $) and (p 2 )^ = ($ I (p'J 2 I $). The latter are converted into the values (r 2 ) and (p 2 ) if 

I $> =| $mtr)- 

Alternatively, according Eqs.(10) we find in the HO model that 



(r 2 )(p 2 ) = (r 2 ) (p 2 )o 



] _ (000 I R I 000) 



_ (000 1 P 1 000) 

ivy; 



(13) 
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Then, taking into account that 

i, 2 3 i, 2 34 

(000 | R I 000) = — r 2 , (000 | P | 000) = — p 2 

and that in the HO model for a nucleus with fully occupied (closed) shells 

(r 2 ) = xr 2 ,, (p 2 ) = xpl, 

x = 2 ^(AT + l)(iV + 2)(iV + 3/2)/A, A = 2 (iV + l)(iV + 2), (14) 



i.e., 



N=0 N=0 



* = § + Z>(^ + + 2 )/ + W N + 2 )» 



AT AT 

where N is the principal quantum number, we get 
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(r\(p\ = x 2 >- A (15) 



<r 2 >(p 2 > = x 2 



A- 1 



2 



1 re - 3/2 

1 + 



1 2 



9 /a - r. 

> 7 — i- (150 



4 V A 



A - I K 

Thus we arrive again to relations (12) and (12'). From this derivation it follows that Eqs. (15) and 
(15') without the signs of inequality are permitted only for the Is 4 configuration in the HO model. 
Such a minimization of the uncertainty relations is retained after making the CM correction. 

3. Further, using the algebraic procedure applied in [5] for the calculation of the expectation 
values (1) and (2) with the Slater determinant |$) = |(ls) 4 ) of the simple HO orbitals, we find 

F(q)=eM-^f)U(q)/U(0), (16) 

U(q) = J exp(-^)/(A,5)dX, 
A 



/(A, q) = (&\ Oi(6) Y[ O a (c) | $) 



a=2 



with the operators 



6 a (x) = exp(— x* ajj exp(— x a a ), (a = 1, . . . , A) 

where b = c + i-^q, c = (^ ~ $ ' ^(^) ^ s the vec tor whose components are the creation 
(annihilation) operators for oscillator quanta in the three different space directions. 
Analogously, one can show that 

S(q, r) = exp(-^r 2 )f/(-^ 2 |r)/f/(0) (17) 



Now, let us assume that the many-body state | $) is a Slater determinant of s.p. states | 
07 ) (7 = 1) • • • ) ^4) which describe completely occupied bound states of the nucleons in a spherically 
symmetric potential well (e.g., Woods-Saxon potential or Hartree-Fock field). Then (see Suppl. B) 

A 

f(X,q)=A- 1 Y,D P , (18) 

where D p is the determinant that is deduced from the determinant 

D=\\ (0 7 I 6(c) I <f) Y ) || 

— * 

replacing the vector c by b in its column with the label p. 

As an illustration, let us consider the (ls) 4 (lp) 12 configuration in the /s-coupling scheme. From 
Eq.(18) it follows that (see Suppl. C) 

f(\,q) = ^d 3 (d 1 + d 2 + d 3 + d i ), (19) 

where d =\\ {4> n im \ 0(c) \ 4>n'i'm') || is the 4x4 determinant, | 4> n i m ) the spatial part of the vector 
I 7 ) in the shell with the radial quantum number n, the orbital angular momentum / and its 
projection m, and di(i = 1, . . . , 4) are deduced from d just as D p from D. 

Subsequent simplifications can be achieved owing to the transformation properties of the matrix 
elements 

M^i m \x) = (<P nlm I O(x) I <t> M ) 
with respect to the rotation group. In fact, we have 

Ml s s {x) = M (x 2 ), Ml; m (x) = M(x 2 )x* m , 

Mg? \x) = M 1 (x 2 )5 mm/ + M 2 (x 2 )x* m x m , (m,m' = -l,0,l) (20) 

with x* = —x. Here x m are the spherical components of the vector x. 

In their turn, the determinants di can be expressed in terms of the scalar functions Mq, M, M\ 
and M2 if one takes into account that the quantities d\ and 0^2 + ^3 + ^4, each separately, are invariant 
under rotations, i.e., they depend on b 2 ,c 2 and be. This property enables us to write down, 

d x = MKc^iM^N^c 2 ) - M(b 2 )M(c 2 )bc}, (21) 

d 2 + d 3 + d 4 = M^iM^Moic 2 ) - M(b 2 )M(c 2 )bc}+ 
2M 1 (b 2 )M 1 (c 2 ){M (c 2 )N (c 2 ) - M 2 (c 2 )c 2 }+ 
+M (c 2 )M 2 (c 2 )M 2 (b 2 )b 2 - M 1 (c 2 )M 2 (b 2 ){M 2 (c 2 ) + M(c 2 )}{b 2 c 2 - (be) 2 } (22) 
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with N {x 2 ) = Afffi 
The scalars M(x 




(23) 



(24) 



from which it follows that 



M,{x 2 ) = exp(r y/4)(A (y) + A 2 (y)), 



(25) 



x 2 M 2 (x 2 ) 



3exp(r V/4)A 2 (|/) 



(26) 



with 



OC 




Jx(yr)Rn(r)r 2 dr, (A = 0, 2) 



o 

where Rn(r) is the radial part of 0n m (r) and j\(z) is the spherical Bessel function of z. In Eqs. 
(25) - (26) x = ir y/V2- 



lation of the simple overlap integrals. 

All these formulae can be useful when studying CM corrections of the cross sections of the elastic 
and quasifree electron scattering from nuclei more complicated than the 4 H e nucleus. An attractive 
feature of similar studies is to proceed with one and the same corrected many-body WF of the 
nuclear g.s. in evaluating different structure functions like the DD and MD. Of certain interest 
might be deviations of the g-dependence of the ratio F(q)/F (q) beyond the HO model from that 
which is given by the canonical Tassie-Barker factor exp(q 2 rl/AA). 
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Appendix A. Comments on simultaneous shrinking of den- 
sity and momentum distributions 

Remind that a commutation relation between two noncommuting operators imposes a definite 
constraint upon their dispersions, i.e., the mean square deviations of these quantities from the 



Thus the initial cumbersome task of handling the expectation value /(A, q) is reduced to calcu 
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corresponding expectation values. In fact, let us consider the two Hermitean operators A and B 
that meet 

A,B] = iC , (A.l) 
where C is also an Hermitean operator. In particular, if A = x and B = p x , then the operator 

c = h 4 . 

By definition, the expectation values with respect to an arbitrary state $ are 



(A) = (ct>L4|ct>>, (B) = ($151$) . 
Let us introduce the deviations 

AA = A-(A), AB = B-(B). 
Obviously they satisfy the relation 



AA,AB 



From Eq.(A.4) it follows the uncertainty relation 



iC . 



1 , - .2 



((AAY)((ABY) > -(C) 



Note also that for operator X, 



((AXf) = ((X-(X)f) = (X 2 )-(X) , 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



if ($ | $) = 1. 

In the case of interest for the relative coordinates r a and the canonically conjugate momenta 
p' a the corresponding relation looks as 

2 



9 A- 1 



4 V A 



or 



<rW> = 9<^><^>> 7 



9 (A- 1 



4V A 

for any state $j nt normalized to unity. Here, in accordance with Eq.(A.6), 

{Q int | (A?) 2 I $mt> = ($mt | (r'i) 2 I = (r 2 ) = J r 2 p(r)dr 



and 



{$ int | (Ap') 2 I *int> = {*int I (Pi)' I *int) = <p") = / P^(p)d P , 



(A. 7a) 



(A.7b) 



(A.8) 



(A.9) 



4 Recall that in our system of units h = 1 
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since (r') = (p') = 0. 

Thus, the general result (A. 5) leads to the condition (A. 7b) for the pair r/ and p^ that obeys 
the commutation rules, 

[(^i)i,(#)J=«W^ (j,* = 1,2,3). (A.10) 
Eqs.(A.lO) are different from 

[(*i)j,(Pi)k]=iS j k U, k — 1,2, 3). (A.ll) 

By the way, this fact means that the transformation r a — > f' a and p Q — > (a = 1, . . . , A) is 
nonunitary, i.e., it cannot be performed via a unitary operator U (UW = IJU^ = 1), viz., putting 
Wr a U = r' a and U^pJJ = p' a . Under the latter the commutation relations would not change. 
It follows from (A. 5) and (A.ll) that 

(r\(p 2 )o > I (A.12) 

for the expectation values 

(r 2 ) = ($ | r\ | $) = y r 2 p (r)dr (A.13) 

(P 2 ) = I Pi I ^> = y P 2 Vo(p)dp (A.14) 

in the shell model ground state $ (in general, any state $) such that ($ | $) = 1. 

So, we have relation (A. 7b) for the mean square radius and mean square momentum in an 
exact translationally invariant state versus relation (A.12) obtained for similar quantities within an 
approximate treatment of the system (say, with the wavefunction which has a violated symmetry 
under space translations). In comparison with Eq.(A.12), the r.h.s. of Eq.(A.7b) is modified by the 
factor {^-) 2 < 1. This modification is closely connected with nonunitarity of the transformation 
from the usual coordinates and canonically conjugate momenta to their relative counterparts. 

Thus the simultaneous shrinking of the density and momentum distributions, shown in the paper 
within the HOM, is consistent with the model independent uncertainty relations (A. 7). 



Appendix B. Comments on derivation of Eq.(18) 

One has to deal with the expectation values of type 

A* = <$|Oi(6)0 2 (c)...6a(c)|^), (B.l) 



and 
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7 (x) = exp(— x*^) exp(xa 7 ) = E^(—x)E 7 (x), (7 = 1, . . . , A) (B.2) 

— * 

where, for instance, b = c + s. One can write 

A* = ($ I E\(-c)El(-c) . . .E^-cjEli-^E^E, (c)E 2 (c) . ..E A (c) | $). (B.3) 
We have used the properties Ei(c+ s) = Ei(c)Ei(s) and [E a (x), Ep(y)\ = for any vectors x and 

y- 

If I $) is a Slater determinant, i.e., 

I $) =| £) e t) = v^4!fi I ... I (B.4) 

with the antisymmetrization operator 

Q = (A\)- 1 J2^ P P (B.5) 
p 

which has the property 

ft 2 = Q, (B.6) 

then 

A$ = (Det(-c) I E^-sjE^s) \ Det(c)} (B.7) 

with 

\Det^) = E 1 (^...E A ^\Det) = ^M 1 (^\M±))---E A (?)\MA)) (B.8) 
Furthermore, using the permutation symmetry of the determinants involved, viz., 

P I Det(-c)) = e p I Det(-c)) (B.9a) 

P I Det(c)) = e p I Det(c)), (B.9b) 

it is easily seen that 

Aq> = (Det(-c) I E^-sjE^s) \ Det(c)} 
= (Det(-c) I E 2 (-s)E 2 (s) \ Det{c)) = ... 

1 A 
= -(Det(-c) I V£ Q (-s)£ a (s) I £>e*(c)) 

Now, taking into account Eq.(£?.6) and the relation 

A 
a=l 
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for any A operators Bi, B 2 , . . . , Ba, we find 

A 



A* = ^i(l) I <$(2) | . . . (MA) | Y.Ei{-Z)E a {s)Y. e v P I ^i(l)> I lfe(2)> • • • | M*)) 



a=l 



where we have introduced the two sets {ip} and {■*/>'} of new orbitals 

| V(a)> = E a (c) | 0(a)) 

and 

| ij;'(a)) = E a (-c) | 0(a)) (a = l,...,A). 
Expression (B.10) explicitly reads 

^ = ^ [ E^^^ 1 ) I ^(-5)^(5) I lMl)><$ I « • • • Wa I V>pJ+ 
p 

+ X> P <# I M)(^ 2 (2) I | lM2)>...<1& I ^p x >+ 

p 

p 

The latter is equivalent to Eq.(18). 



(B.10) 
(B.lla) 
(B.llb) 



Appendix C. Evaluation of determinants in Eq.(18) 

In order to simplify evaluation of the (ls) 4 (lp) 12 configuration determinants involved in the r.h.s. 
of Eq.(18), let us consider a sparse nm x nm matrix 

Zu(m x to) Z 12 (m x to) ... Z ln (m x to) 
Z 21 (m x to) Z 22 (m x to) ... ^2n( m x m ) 



Z{nm x nm) 



Z n i(m x to) ^„2( m x to) ... Z nri (m x to) 



that consists of n 2 diagonal mxm block matrices 



Z tt (roxm)=diag[Z<^>,...,Z< 



'Hi 
ik J 



or 



Z ik (m x to) 



• ... 
• ... 

... • 



(C.l) 



(C.2a) 



(C.2b) 
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(i,k = l,2,...,n). 

In the representation (C.2b) the diagonal elements zf^ (1=1,2,. . . ,m) which are, in general, arbitrary 
and different from one another, are marked by the symbol •. One can show that det Z =|| Z{nm x 
nm) || is equal to the product 

m 

det Z = Y[detYi (C.3) 
i=i 

of the m determinants of the nxn matrices Yj, (i = 1, . . . ,m). Each of Yi is composed of the n 2 
diagonal elements of the matrices Z ik (m x m) (namely, elements •), which one encounters when 
passing in the clockwise (or counter-clockwise) direction the corresponding spiral-like contour as 
displayed in Fig. 2. 

One should note that these m contours are not closed (i.e., they are open), and the i-th contour 
begins at the element z[^ on the diagonal of the block Z n . 

The relationship (C.3) can be proved using the simple properties of determinants. As an illus- 
tration, let us demonstrate relation (C.3) for m = 2 and n = 3: 
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Figure 2: The relevant contours in the case with m = 2 and n = 3. 



Therefore 

ai bi 
a 2 
Ci d\ 
c 2 
5-i /ii 
g 2 

The aforementioned prescription yields the same result but much quicker. 
In the case of interest one has to deal with the matrices of the type 



M = 



M ss 


M s+ 


M s0 


M s _ 


M +s 


M ++ 


M +0 


M+_ 


M 0s 


M Q+ 


M 00 


M _ 




M„+ 


M_ 


M__ 



where each of the sixteen 4x4 blocks is diagonal, e.g., 

M ss = diag[Mi, Ml, Ml Mi] 
with the four equivalent dispositions: 

diag[M£$, M^(c1, M£@, M^c)], 

diagfM^c), M^(6), Mi'(g), M£(a)], 

diagfM^c), M£(S), M£(S), M£(<9], 

diag[M£(S), Ml:(c), Ml:(c), M%@)]. 

Taking into account this equivalence we get expression (19) which is much simpler compared to 
Eq.(18). 
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b 2 e 2 

ft 

d 2 h 
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/l2 j 2 



ai 6i ei 
ci di /i 
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c 2 d 2 / 2 
#2 /la J2 
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